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The directed flow of protons and 7r + have been studied in 158 A GeV 208 Pb + 208 Pb collisions. 
A directed flow analysis of the rapidity dependence of the average transverse momentum projected 
onto the reaction plane is presented for semi-central collisions with impact parameters m 8 fm, where 
the flow effect is largest. The magnitude of the directed flow is found to be significantly smaller 
than observed at AGS energies and than RQMD model predictions. 

25.75.-q,25.75.Ld,25.75.Dw 



1 



Collective flow has been studied in heavy ion colli- 
sions since first observed at the Bevalac by the Plastic 
Ball experiment Q. At Bevalac energies of a few GeV 
per nucleon and lower, the study of collective flow has 
been of interest largely due to its expected sensitivity 
to the nuclear equation of state (EOS) However, 
the extraction of information on the EOS in heavy ion 
collisions is complicated by uncertainties in the initial 
dynamics of the pre-hydrodynamic stage, such as due 
to the momentum-dependence of the repulsive nucleon- 
nucleus interaction |3j and possible in-medium modifi- 
cations of the nucleon-nucleon cross section 0. Never- 
theless, the importance of collective flow measurements 
in ultra-relativistic heavy ion collisions has been empha- 
sized by several authors Collective flow develop- 
ment follows the time evolution of pressure gradients in 
the hot, dense matter. Thus, collective flow can serve 
as a hadronic "penetrating probe" to provide informa- 
tion on the initial state. In particular, the formation of 
a Quark Gluon Plasma (QGP) during the early stages 
of the collision is expected to result in reduced pressure 
gradients due to a softening of the EOS, with a corre- 
sponding reduction of collective flow • 

Transverse collective flow is normally discussed in 
terms of its lowest order symmetries with respect to the 
reaction plane, which have recently been formulated in 
terms of a Fourier decomposition jl(|[ll| ■ The lowest or- 
der component is radial flow which is characterized by 
an isotropic transverse flow velocity. The next compo- 
nent is the directed flow which is characterized by the 
net displacement of the flow into a particular transverse 
direction. The elliptic flow component corresponds to the 
second order Fourier coefficient of the flow pattern. 

Recently, collective flow has been observed at the AGS 
energy of 11 A GeV 0|. The directed flow is observed 
to be smaller than at lower incident energy but similarly 
consistent with model calculations jyj). It is also ob- 
served that the elliptic flow has changed from an out-of- 
plane squeeze-out direction to an in-plane direction [fl3| . 
At CERN SPS energies it has recently been shown that 
it is possible to determine the event-plane in the tar- 
get fragmentation region for 16 O- and 32 S-induced reac- 
tions jT|] and also in the mid-rapidity region for 32 S- fT^ j 
and 208 Pb-induced reactions jl6| . Directed flow results at 
SPS energies were first reported in |l7| and both directed 
and elliptic flow results have recently been published jl8| . 
In this letter we analyze the centrality dependence of the 
directed flow of protons and ir + in 158 A GeV 208 Pb + 
208 Pb and present a detailed analysis of its rapidity de- 
pendence at the intermediate centrality where it is great- 
est. 

The present analysis makes use of a subset of the de- 
tector systems of the WA98 experiment. This subset 
consists of the trigger detectors, the Plastic Ball detec- 
tor, and the tracking spectrometers. The large aperture 
dipole magnet, Goliath, provided momentum analysis for 



the tracking detectors. The data presented were taken 
during the 1996 SPS run period with 158 A GeV 208 Pb 
beams using a 213 /jm thick 208 Pb target. The WA98 
minimum bias cross section for this run period, with mag- 
netic field on, was a m b = 6450 mb. 

The trigger detectors consisted of a nitrogen gas 
Cerenkov counter to provide a fast beam trigger (< 30 ps 
time resolution), beam- halo veto counters, and the 
MIR AC calorimeter. A beam trigger was defined as a 
signal in the start counter with no coincident signal in the 
veto counter (which had a 3 mm diameter circular hole) 
or in beam halo counters. The MIRAC measures the to- 
tal transverse energy over the interval 3.2 < r\ < 6.0 with 
full azimuthal coverage over the interval 3.7 < rj < 4.9. 
The WA98 minimum bias trigger requires a clean beam 
trigger with a MIRAC transverse energy which exceeds 
a low threshold. 

The Plastic Ball detector provides full azimuthal cover- 
age in the target fragment region (pseudorapidity — 1.7 < 
rj < 0.5) with 655 detector modules. It provides iden- 
tification of pions, protons, deuterons, and tritons (tt, 
p, d, and t) with kinetic energies of 50 to 250 MeV 
by the AE — E method. In addition, stopped ir + are 
identified by detection of the delayed e + from the de- 
cay 7r + — > /i + + — > e + + v e + Vfj_ + Vfj, in the Plas- 
tic Ball. For the present analysis the rapidity region 
—0.6 < y(proton) < 0.3 has been used. 

The measurement of identified particles near mid- 
rapidity was obtained using two large acceptance track- 
ing arms beginning about 3.3 meters downstream of the 
Goliath magnet. With the normal operation magnetic 
field setting the first tracking arm was located to the 
side of negative charge deflection and the second tracking 
arm to the side of positive charge deflection. The posi- 
tive charge results are presented here. The momentum 
resolution of the tracking system may be parameterized 
as Ap/p ~ 0.97% + 0.16% p + 0.023% p 2 {p in GeV/c). 
The acceptance for protons covered an interval Ay p 0.3 
which shifted with transverse momentum to provide cov- 
erage over the region IA < y p < 2.4. Particle identifica- 
tion was obtained by time-of-flight measurement with a 
resolution of < 90 ps. 

We determine the reaction plane as the azimuthal di- 
rection, $o, opposite to Pt, the total transverse momen- 
tum vector of fragments (p, d, and t) detected in the 
target rapidity region in the Plastic Ball detector.[] To 
check for detector effects, mixed events are created by 
mixing particles from different events keeping the same 
multiplicity distribution as in the real events. The mixed 
events are then analyzed in the same manner as the real 
events. The laboratory distribution of <i>o is nearly uni- 



1 Note: the opposite direction is chosen by convention that 
the projectile fragments define the direction of positive flow. 



form, with less than 2% variations. The $ mtx distribu- 
tion for mixed events is observed to show the same weak 
variations indicating small detector effects. The data 
have not been corrected for this effect since the analy- 
sis of the mixed events shows a negligible result. 

In order to study how well the fragment flow direc- 
tion is defined, we divide each event randomly into two 
equal sized subevents and determine a fragment flow di- 
rection for each subevent, <3? a and If the direction 
of Pt is well-defined, the directions determined from 
each subevent should be strongly correlated [p pS|JTl[ | . In 
Fig. |l|, the <I> a — (f>b correlation is shown for two different 
centrality bins. As expected, the correlation observed for 
semi-central events is significantly larger than for very 
central events. Also shown are the results for mixed 
events allowing (open squares) or forbidding (open cir- 
cles) multiple module hits. The former case demonstrates 
that detector non-uniformities are negligible. In the lat- 
ter case a weak anti-correlation is observed due to the 
finite detector granularity and an excluded-module ef- 
fect. The data have been corrected for this effect using 
such mixed events. 
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1. The distribution of differences between the to- 



tal transverse momentum directions of two randomly chosen 
equal size subevents of fragments (p,d,t) in the target rapidity 
region for a) semi-central (100 < Et < 200 GeV) and b) cen- 
tral (380 <E T < 420 GeV) collisions of 158 A GeV 208 Pb + 
208 Pb. Solid circles are for subevents within the same event. 
Open points are for subevents constructed from mixed events. 
Solid curves are fits to guide the eye. 



Azimuthal anisotropics of the particle emission are 
evaluated by means of a Fourier expansion pO|p] |. The 
Fourier coefficients v n (n = 1,2) are extracted from the 
azimuthal distribution of identified particles with respect 
to the reaction plane, which is determined using all 
other fragments in the Plastic Ball. 
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where ef> is the measured azimuthal angle. The Fourier co- 
efficient v[ quantifies the directed flow, whereas v' 2 quan- 



tifies the elliptic flow. The coefficients must be corrected 
for the event plane resolution as v n = v' n / (cos(n(&o — 
<F r ))} where $o — &r is the deviation of the measured 
reaction plane from the true reaction plane. The event 
plane resolution may be extracted from the correlation 
between subevents. For weak correlations one expects 
(cos($ - $r)) — v/2( cos ( < J > o - $&))• Using the more ac- 
curate procedure and interpolation formula of Ref. p"l| ] 
one obtains (cos($o - $r)) = 0.377 ±0.018 for the semi- 
central (100 < Et < 200 GeV) event selection. 

The dependence of the vi fit parameter on centrality, 
as determined by the measured transverse energy (Et), 
is shown in Fig. ||. For convenience an impact param- 
eter scale is also shown. The Et scale has been con- 
verted to an impact parameter scale assuming a mono- 
tonic relationship between the two quantities, and equat- 
ing da/dET with da/db. As seen in Fig. ||, the strength of 
the directed flow of protons increases with centrality and 
reaches a maximum value for semi-central collisions with 
b w 8 fm. It is interesting to note that the strongest flow 
effect occurs at larger impact parameters than observed 
at lower incident energy for similar systems (where b w 4 
fm) 0,0. For comparison, RQMD 2.3 @ model pre- 
dictions are shown subjected to the same analysis after 
applying the Plastic Ball detector acceptance, but using 
the true reaction plane. RQMD predicts a significantly 
stronger correlation for protons than observed. 
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FIG. 2. The centrality dependence of the directed flow co- 
efficient Vi for protons (circles) and ir + (squares). Triangles 
are results from RQMD model calculations. The data have 
been corrected for the event-plane resolution. The vertical 
^ jaars indicate the uncertainty of the fit and resolution correc- 
tion. The horizontal bars indicate the Et bin intervals (or 
impact parameter intervals for RQMD). 



Also shown in Fig. [2] is the strength of the directed 
flow of 7r + , identified in the Plastic Ball. A clear anti- 
correlation, or anti-flow ^0|, is observed between the 
fragment and 7r + flow directions. This behaviour has 
been observed at incident energies from 1 A GeV to SPS 
energies and has been explained as resulting from prefer- 
ential absorption of the pions emitted in the target spec- 
tator direction [ |l^ , |l4| , p0| , pl| . The absorption results in 
an oppositely directed apparent n + flow. The strength 
of the anti-correlation increases for the most peripheral 
events, indicating the increasing role of absorption. 

A conventional directed flow analysis has been per- 
formed Jl9[| , in which the average transverse momentum 
with respect to the reaction plane (p x ) is evaluated as a 
function of rapidity. This is done for semi-central colli- 
sions (100 < Et < 200 GeV) where the largest azimuthal 
asymmetry is observed (see Fig. ||). The distribution 
d 3 N / dp x dp' y dy is constructed for protons and ir + in the 
Plastic Ball and in the tracking arm, where the new axis 
p' x corresponds to the reaction plane determined event- 
by-event using all remaining fragments measured in the 
Plastic Ball (then reflected, p' x — > —p' x , to correspond to 
the projectile fragment direction, according to conven- 
tion). At each rapidity the average transverse momen- 
tum in the reaction plane, {p' x ), is calculated from fits to 
the experimental distributions. 

Similar to v n , the average projected momenta are re- 
duced by (p' x ) = (p x ) ■ (cos($o — $»■))< After correction 
for the event-plane resolution, the (p x ) for protons and 
7r + are plotted as a function of rapidity in Fig. |^. As 
expected from Fig. pi the n + show an anti-flow relative 
to the proton flow. Summing over the Plastic Ball accep- 
tance (p x ) values of 8.2 ± 0.7, -24.9 ± 1.9, -53.6 ± 4.1, 
and —78.4 ± 5.8 MeV/c are obtained for 7r + , p, d, and 
t, respectively. The observed scaling with fragment mass 
for p, d, and t indicates emission sources with a common 
collective motion. 

The main sources of systematic error in the present 
analysis are: detector non-uniformities, contamination 
in the particle identification, and fit biases in extract- 
ing (p x ) [p3| . The mixed event analysis indicates that 
systematic errors from detector non-uniformities are less 
than 2%. The effect of contamination has been estimated 
by Monte Carlo simulations. The amount of contamina- 
tion can be estimated from fitting the background under- 
lying the peaks in the Plastic Ball particle identification 
spectra. For example, the amount of contamination in 
the proton sample varies from 6% in peripheral events to 
26% in central events. The effect of contamination has 
been estimated in simulation by analyzing events with 
various amounts and types of contaminated particle dis- 
tributions compared to pure particle distributions, where 
the various distributions are taken to have spectra and 
flow characteristics similar to those measured. These 
studies indicate a maximum systematic error of 8.5%. 
The extraction of (p x ) is estimated to have an additional 



15% uncertainty deduced from observed variations in the 
results depending on the fit region or method used to fit 
the d 3 N I dp' x dp' y dy distribution. These systematic errors 
have not been included in Figs. ^ or|3[ 
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FIG. 3. The average transverse momentum projected onto 
the reaction-plane for semi-central 158 A GeV Pb + Pb 
collisions (note y cm = 2.9). The vertical errors indicate the 
statistical errors of the fit only. The horizontal bars on the 
tracking points indicate the width of the rapidity bin. RQMD 
model calculations (b=8-10 fm) and VENUS model calcula- 
tions (b=8-10 fm) are also shown. The VENUS prediction for 
7r + (not shown) is similar to that of RQMD. 



In Fig. |3J the measured results are compared to 
RQMD 2.3 H and VENUS 4.12 || predictions for 
similar impact parameter range. The RQMD calcula- 
tion, in cascade mode, overpredicts the observed proton 
flow by about a factor of three. On the other hand, at 
AGS energies cascade mode RQMD calculations under- 
predict the observed directed flow by about a factor of 
two, but reasonable agreement is obtained when mean 
field effects are included Q|. At SPS energies mean field 
effects are expected to be smaller, but would increase the 
observed disagreement. The VENUS predictions show a 
similar disagreement in the target rapidity region. The 
results suggest a significant softness in the nuclear re- 
sponse. The maximum proton (p x ) observed is in better 
agreement with predictions of the Quark Gluon String 
Model (with rescattering) of Ref. | ] and with a 3-fluid 
hydrodynamical model calculation [25|. However, these 
predictions have not been filtered with the experimental 
acceptance and both calculations predict that the max- 
imum (p x ) occurs about one unit forward of the target 
rapidity. It is interesting to note that VENUS predicts 
a complicated proton flow behaviour with protons hav- 



ing an anti-flow direction (similar to the RQMD pion 
prediction) near mid-rapidity. However, this prediction 
disagrees with the results of Ref . |Q . 

In summary, the directed flow of protons and 7r + has 
been studied in 158 A GeV 208 Pb + 208 Pb collisions. 
The directed flow is largest for impact parameter rs 8 
fm, which is considerably more peripheral than observed 
at lower incident energies. The tt + directed flow is in the 
direction opposite to the protons, similar to observations 
at 11 A GeV energy The magnitude of the pro- 

ton directed flow is much less than cascade mode RQMD 
model predictions, which underpredict the proton flow at 
AGS energies. It is also much less than VENUS model 
predictions. The results indicate a soft nuclear response 
compared to these model predictions at SPS energies. 

We express our gratitude to the CERN accelerator di- 
vision for the excellent performance of the SPS accelera- 
tor complex. We gratefully acknowledge the effort of all 
engineers, technicians, and support staff who have made 
possible the construction and operation of this experi- 
ment. 

This work is supported by the German BMBF and 
DFG, the U.S. DOE, the Swedish NFR and FRN, 
the Dutch Stichting FOM, the Department of Atomic 
Energy, the Department of Science and Technology, 
and the University Grants Commission of the Govern- 
ment of India, the Indo-FRG Exchange Programme, 
the Stiftung fur Deutsch-Polnischc Zusammenarbeit, the 
Grant Agency of the Czech Republic under contract No. 
202/95/0217, the PPE division of CERN, the Swiss Na- 
tional Fund, the INTAS under contract INTAS-97-0158, 
the Grant-in- Aid for Scientific Research (Specially Pro- 
moted Research & International Scientific Research) of 
the Ministry of Education, Science and Culture, JSPS 
Research Fellowships for Young Scientists and also by 
the University of Tsukuba Special Research Projects, 
and ORISE. ORNL is managed by Lockheed Martin 
Energy Research Corporation under contract DE-AC05- 
960R22462 with the U.S. Department of Energy. MIT is 
supported by the U.S. Department of Energy under the 
cooperative agreement DE-FC02-94ER40818. 



[9 
[10 

[11 

[12; 

[13; 
[14; 
[15; 

[16 

[17; 

[is: 

[19 

[20 
[21 
[22 
[23 

[24 
[25 



(1995). 

D.H. Rischke, Nucl. Phys. A610 (1996) 88c. 

H. Sorge, Phys. Rev. Lett. 78, 2309 (1997). 

S.A. Voloshin and Y. Zhang, Z. Phys. C70, 665 (1996). 

A.M. Poskanzer and S. A.Voloshin, Phys. Rev. C58, 1671 

(1998). 

E877 Collaboration, J. Barrette, et al, Phys. Rev. Lett. 
73, 2532 (1994). 

E877 Collaboration, J. Barrette, et al, Phys. Rev. C 56, 
3254 (1997). 

WA80 Collaboration, T.C. Awes, et al, Phys. Lett. B 
381, 29 (1996). 

WA93 Collaboration, M.M. Aggarwal, et al, Phys. 
Lett. B 403, 390 (1997). 

NA49 Collaboration, T. Wienold, et al, Nucl. Phys. A 
610, 76c (1996). 

WA98 Collaboration, M. Kurata, et al, p. 549 and 
S. Nishimura, et al, p. 258, Physics and Astrophysics of 
Quark Gluon Plasma, Narosa Publishing House (1998). 
NA49 Collaboration, H. Appelshauser, et al, Phys. Rev. 
Lett. 80, 4136 (1998). 

P. Danielewicz and G. Odyniec, Phys. Lett. 157B, 146 
(1985). 

A. Jahns, et al, Phys. Rev. Lett. 72, 3463 (1994). 

A. Kugler, et al, Acta Phys. Pol. 25, 691 (1994). 

H. Sorge, Phys. Rev. C 52, 3291 (1995). 

M. Kurata, PhD thesis, University of Tsukuba. 

H. Schlagheck, PhD thesis, University of Miinster. 

K. Werner, Phys. Rep. 232, 87 (1993). 

A. Dumitru, et al, Heavy Ion Phys. 5, 357 (1997). 



[1] H.A. Gustafsson, et al, Phys. Rev. Lett. 52, 1590 (1984). 
[2] H. Stocker and W. Greiner, Phys. Rep. 137, 277 (1986). 
[3] J. Aichelin, et al, Phys. Rev. Lett. 58, 1926 (1987). 
[4] See for example, W. Reisdorf and H.G. Ritter, Ann. Rev. 

Nucl. Part. Sci. 47, 663 (1997). 
[5] N.S. Amelin, et al, Phys. Rev. Lett. 67, 1523 (1991). 
[6] J.-Y. Ollitrault, Phys. Rev. D 46, 229 (1992); J.-Y. Olli- 

trault, Phys. Rev. D 48, 1132 (1993). 
[7] CM. Hung and E.V. Shuryak, Phys. Rev. Lett. 75, 4003 



